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We reanalyze deep inelastic scattering data of BCDMS Collaboration by including proper cuts 
of ranges with large systematic errors. We perform also the fits of high statistic deep inelastic 
scattering data of BCDMS, SLAC, NM and BFP Collaborations taking the data separately and 
in combined way and find good agreement between these analyses. We extract the values of both 
the QCD coupling constant a s (M§) up to NLO level and of the power corrections to the structure 
J> ' function F%- 

oo : 
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■ 1. Introduction 

00 ' 

The deep inelastic scattering (DIS) leptons on hadrons is the basical process to study the 
CN| ■ values of the parton distribution functions (PDF) which are universal (after choosing of fac- 

torization and renormalization schemes) and can be used in other processes. The accuracy of 
the present data for deep inelastic structure functions (SF) reached the level at which the Q 2 - 
dependence of logarithmic QCD-motivated terms and power-like ones may be studied separately 
O j' (for a review, see [1] and references therein). 

In the present paper we sketch the results of our analysis [2] at the next-to-leading order 
(NLO) of perturbative QCD for the most known DIS SF F*2(x,Q 2 ) taking into account experi- 
mental data [4]- [7] of SLAC, NM, BCDMS and BFP Collaborations. We stress the power-like 
effects, so-called twist-4 (i.e. ~ l/Q 2 ) contributions. To our purposes we represent the SF 
F 2 (x, Q 2 ) as the contribution of the leading twist part F% (x, Q 2 ) described by perturbative 
QCD, when the target mass corrections are taken into account, and the nonperturbative part 
("dynamical" twist-four terms): 

F 2 (x, Q 2 ) = Ft\x, Q 2 ) = Ff CD {x, Q 2 ) (l + , (1) 

where h&(x) is magnitude of twist-four terms. 

Contrary to standard fits (see, for example, [8]- [10]) when the direct numerical calculations 
based on Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equation [11] are used to evaluate 
structure functions, we use the exact solution of DGLAP equation for the Mellin moments 
Mn w2 (Q 2 ) of SF F 2 tw2 (x,Q 2 ) and the subsequent reproduction of F 2 /u "(x,Q 2 ) at every needed 
Q 2 -v&lue with help of the Jacobi Polynomial expansion method 1 [12, 13] (see similar analyses 
at the NLO level [13, 14] and at the next-next-to-leading order (NNLO) level and above [15]. 

In this paper we do not present exact formulae of Q 2 -dependence of SF F 2 which are given 
in [2]. We note only that the PDF at some Qq is theoretical input of our analysis and the 
twist-four term h&{x) is considered as a set of free parameters (one constant h^ixi) per Xj-bin): 
h^ ee {x) = J2i=i h±{xi) } where I is the number of bins. 



* Talk presented at X Int. Workshop on Deep Inelastic Scattering (DIS2002), Cracow, 30 April - 4 May 2002. 

A shorter version will appear in the Proceedings. 
1 We note here that there is similar method [16], based on Bernstein polynomials. The method has been used 

in the analyses at the NLO level in [17] and at the NNLO level in [18]. 
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2. Fits of F 2 

First of all, we choose the cut Q 2 > 1 GeV 2 in all our studies. For Q 2 < 1 GeV 2 , the 
applicability of twist expansion is very questionable. Secondly, we choose Qq = 90 GeV 2 (Qq = 
20 GeV 2 ) for the nonsinglet (combine nonsinglet and singlet) evolution, i.e. quite large values 
of the normalization point Qq: our perturbative formulae should be applicable at the value of 
Qq. Moreover, the higher order corrections ~ cc^Q 2 ,) and ~ (a s (Q 2 ) — a s (Q'o)) k (k > 2) should 
be less important at these Qq values. 

We use MINUIT program [19] for minimization of \ 2 {F2). We consider free normalizations of 
data for different experiments. For the reference, we use the most stable deuterium BCDMS data 
at the value of energy Eq = 200 GeV (Eq is the initial energy lepton beam). Using other types 
of data as reference gives negligible changes in our results. The usage of fixed normalization for 
all data leads to fits with a bit worser x 2 - 

2.1. BCDMS 12 C + H 2 + D 2 data 

We start our analysis with the most precise experimental data [6] obtained by BCDMS muon 
scattering experiment at the high Q 2 values. The full set of data is 762 points. 

It is well known that the original analyses given by BCDMS Collaboration itself (see also 
Ref. [9]) lead to quite small values a s (M§) = 0.113. Although in some recent papers (see, for 
example, [8, 20]) more higher values of the coupling constant q s (M|) have been observed, we 
think that an additional reanalysis of BCDMS data should be very useful. 

Based on study [21] we proposed in [2] that the reason for small values of a s (M|) coming 
from BCDMS data was the existence of the subset of the data having large systematic errors. 
We studied this subject by introducing several so-called V-cuts 2 (see [2]): 

y > 0.14 when 0.3 < x < 0.4, y > 0.16 when 0.4 < x < 0.5 
y > V cut 3 when 0.5 < x < 0.6, y > Y cut 4 when 0.6 < x < 0.7 
y > Y cut5 when 0.7 < x < 0.8 (2) 

and several N sets for the cuts at 0.5 < x < 0.8: 



N 





1 


2 


3 


4 


5 


6 


Y C ut3 





0.14 


0.16 


0.16 


0.18 


0.22 


0.23 


Y C utA 





0.16 


0.18 


0.20 


0.20 


0.23 


0.24 


Y C ut5 





0.20 


0.20 


0.22 


0.22 


0.24 


0.25 



Table 1. The values of Y cut3 , Y cut4 and Y cut5 . 



From the Figs. 1 and 2 we can see that the a s values are obtained for iV = 1 6 of Y cut 3, 
YcutA and Y cut ^ are very stable and statistically consistent. 

2.2. SLAC, BCDMS, NM and BFP data 

After these V-cuts have been incorporated (with = 6) for BCDMS data, the full set of 
combine data is 1309 points. The results of the fits are compiled in Summary. 

3. Summary 

We have demonstrated several steps of our study [2] of the Q 2 -evolution of DIS structure 
function F 2 fitting all fixed target experimental data. 

From the fits we have obtained the value of the normalization a s (M§) of QCD coupling con- 
stant. First of all, we have reanalyzed the BCDMS data cutting the range with large systematic 
errors. As it is possible to see in the Figs. 1 and 2, the value of a s (M^) rises strongly when the 

2 Hereafter we use the kinematical variable Y = (Eo — E)/Eo, where E is scattering energies of lepton. 
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Fig. 1. The study of systematics at different Y cut 
values in the fits based on nonsinglet evolution (i.e. 
when x > 0.25). The inner (outer) error-bars show 
statistical (systematic) errors. 



Fig. 2. All other notes are as in Fig. f with two 
exceptions: the fits based on combine evolution 
and the points Ny C ut = 1, 2, 3, 4, 5 correspond the 
values N =1,2, 4, 5, 6 in the Table 1. 
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Fig. 3. The values of the twist-four terms. The black and white points correspond to the small- a; asymp- 
totics ~ of sea quark and gluon distributions with w = and w = 0.18, respectively. The statistical 
errors are displayed only. 



cuts of systematics were incorporated. In another side, the value of a s (M|) does not dependent 
on the concrete type of the cut within modern statistical errors. 



4 



DIS02L PRINTED ON FEBRUARY 7, 2008 



Fitting SLAC, BCDMS, NM and BFP data, we have found in [2] that at Q 2 > 10-=- 15 GeV 2 
the formulae of pure perturbative QCD (i.e. twist-two approximation together with target mass 
corrections) are in good agreement with all data 3 . When we have added twist-four corrections, 
we have very good agreement between QCD (i.e. first two coefficients of Wilson expansion) and 
the data starting already with Q 2 = 1 GeV 2 , where the Wilson expansion should begin to be 
applicable. The results for a s (M§) are very similar (see [2]) for the both types of analyses and 
have the following form: 

a s {M 2 z ) = 0.1177 ±0.0007 (stat) ±0.0021 (syst) ± 0.0009 (norm), (3) 

where the symbols "stat" , "syst" and "norm" mark the statistical error, systematic one and the 
error of normalization of experimental data. 

We would like to note that we have good agreement also with the analysis [20] of combined 
HI and BCDMS data, which has been given by HI Collaboration very recently. Our results for 
a s (M§) are in good agreement also with the average value for coupling constant, presented in 
the recent studies (see [8, 28] and references therein) and in famous reviews [29]. 

At the end of our paper we would like to discuss the contributions of higher twist corrections. 
In our study [2] we have reproduced well-known x-shape of the twist-four corrections at the large 
and intermediate values of Bjorken variable x (see the Fig. 3 and [9, 2]). 

Note that there is a small- x rise of the magnitude of twist-four corrections, when we use flat 
parton distributions at x — > 0. As we have discussed in Ref [2], there is a strong correlation 4 
between the small- x behavior of twist-four corrections and the type of the corresponding asymp- 
totics of the leading-twist parton distributions. The possibility to have a singular type of the 
asymptotics leads (in our fits) to the appearance of the rise of sea quark and gluon distributions 
as ~ x -°- 18 a t low x values, that is in full agreement with low x HERA data and with theoretical 
studies [22, 26]. At this case the rise of the magnitude of twist-four corrections is completely 
canceled. This cancellation is in full agreement with theoretical and phenomenological studies 
(see [22, 31, 23]). 
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